Can Kinematic Diffraction Distinguish Order from Disorder? 
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Diffraction metfiods are at tiie iieart of structure determination of solids. Wfiile Bragg-Iike scat- 
tering (pure point diffraction) is a cliaracteristic feature of crystals and quasicrystals, it is not 
straightforward to interpret continuous diffraction intensities, which are generally linked to the 
presence of disorder. However, based on simple model systems, we demonstrate that it may be 
impossible to draw conclusions on the degree of order in the system from its diffraction image. In 
particular, we construct a family of one-dimensional binary systems which cover the entire entropy 
range but still share the same purely diffuse diffraction spectrum. 

PACS numbers: 61.05.cc 61.43.-j 61.44.-n 02.50.-r 
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The inverse problem of reconstructing a structure from 
its diffraction pattern is one of the most important chal- 
lenges in materials science. Its degree of complexity in- 
creases if one goes beyond simple periodic systems to 
cover quasicrystals, modulated structures or complex al- 
loys. In particular, it has been realised that the Bragg 
diffraction alone is generally insufficient for a complete 
reconstructionii^. 

Currently, an increasing effort is being made to under- 
stand and utilise the continuous part of the diffraction 
pattern; see Refs. for background and Refs. UiSj for 
recent applications. However, even in the idealised situ- 
ation of a perfect diffraction experiment with unlimited 
resolution, the reconstruction is generally not unique. Al- 
ready in 1944, Patterson^ discussed homometric point 
sets, which are point sets whose kinematic diffraction 
patterns coincide, and provided explicit examples to il- 
lustrate the ambiguity. It was demonstrated^ that it may 
be possible to lift the ambiguity, and thus to determine 
the structure uniquely, if higher-order correlations are 
known. While one can argue that, for structures origi- 
nating from systems with pure pair-potential interaction 
(or allowing a description by effective pair potentials, 
compare Ref.0), higher-order correlations are determined 
by the pair correlations^, this is not generally the case, 
and in practice measurements of higher-order correlations 
are extremely difficult. The role of phase information in 
stochastic systems was investigated in Ref. 9. 

Here, we want to go one step further, and compare the 
diffraction of various point sets, ranging from determinis- 
tic to fully stochastic, in a parametrised way. We charac- 
terise the degree of order by the corresponding (metric) 
entropy. As we will see below, it is possible to construct 
families of point sets which span an entire entropy range 
but share the same kinematic diffraction pattern - prov- 
ing that diffraction is insensitive even to the degree of 
order in this case. 

We start by giving a brief introduction to some ba- 
sic notions of mathematical diffraction theory, for a one- 



dimensional (but relevant) setting with scatterers placed 
on integer positions. Although this is a highly idealised 
situation that ignores displacement effects, its practi- 
cal relevance is well known; see Ref. [l| and references 
therein. Then, we discuss the diffraction of two binary 
systems (characterised by two scattering strengths) - a 
perfectly ordered structure based on a specific determin- 
istic sequence, and a completely random structure based 
on a coin-tossing experiment. It was observed earlier— 
that these rather different systems share the same diffrac- 
tion, and are thus homometric. Finally, we introduce a 
'Bernoullisation' procedure to couple coin-tossing disor- 
der to a perfectly ordered structure, thus producing par- 
tially ordered systems of varying entropy. We employ this 
procedure to explicitly construct a family of binary sys- 
tems which are homometric and cover the entire available 
entropy range from the perfectly ordered (entropy 0) to 
the fully stochastic situation (entropy log(2)). Although 
these systems may not occur naturally, they can be made 
synthetically. In the simplest scenario, a binary struc- 
ture can be produced by sequential deposition of layers 
consisting of two different materials, which makes it pos- 
sible to realise any desired sequence; see Ref. [ll| for an 
example. More complicated structures are also feasible, 
see for instance Ref. 12, and such artificial materials with 
designed physical properties will become increasingly im- 
portant. 

a. Diffraction of Dirac combs. To keep arguments 
simple, we consider the diffraction of one-dimensional 
systems with point-like scatterers located at integer 
points n e Z. The scattering strengths are given by 
weights Wn for n G Z, which we assume to be real for sim- 
plicity (the setting can be extended to complex weights). 
The corresponding scattering density is modelled by the 
Dirac comb 
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where 5x denotes the normahsed point measure (Dirac 
5) on the real line, located at position x. Clearly, all 
distances between scatterers are integer valued. This im- 
plies that the autocorrelation (or Patterson) measure 7, 
assuming its existence for the moment, is again a Dirac 
comb on Z, 
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with the coefficients 77(771) obtained as the limits 

1 ^ 

rj(m) — lim > WnWn+m- 

The scattering intensity I{k) for wave numbers fc e M is 
then determined by the diffraction measure 7, the Fourier 
transform of the autocorrelation 7, compare Ref. [l^ for 
background. There are several slightly different versions 
of the Fourier transform. We prefer to use 
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for a Schwartz function 0, and its standard extension 
to tempered distributions and measures; see Ref. [ij for 
details. 

For the case of a one-dimensional crystal with Wn = w 
for all 77, G Z, we have 77(777) = uP' for all 777 G Z, hence 
7 = vj^ 5j^, where we use (5^ as shorthand for the sum 
X^nGZ'^"- Fourier transform is obtained by Poisson's 
summation formula^, & = which gives 7 = 6z- 
The diffraction image thus consists entirely of Bragg 
peaks, located at integer positions k, with equal diffrac- 
tion intensities I{k) = . The diffraction spectrum in 
this case is pure point, meaning that it consists of Bragg 
peaks only. In general, the diffraction measure may com- 
prise three different contributions, 

7 = 7pp + 7sc + 7ac 

where 7pp is the pure point part, consisting of a countable 
sum of 5 peaks. The term 7ac corresponds to the abso- 
lutely continuous component, which can be described by 
a locally integrable (and often continuous) non-negative 
function /ac(fc) of the wave vector k. The remainder, if 
there is any, is called the singular continuous component 
75C. While it vanishes on the complement of a set S of 
measure 0, even within S it never gives weight to any 
single point. When such a component is present, S can 
thus not be a countable set. Apart from trivial examples 
of a diffraction measure that is concentrated on a line in 
the plane, or similarly on a manifold of lower dimension, 
typical examples for this strange contribution are diffrac- 
tion intensities which are supported on a Cantor set or 
a dense set. A well-known example for the latter phe- 
nomenon is the Thue-Morse chain. For an appropriate 
choice of its scattering strengths, it has a purely singu- 
lar continuous diffraction measure; see Refs. 1113171 for 
derivations and Refs. [TllfT3 for applications. 



If the Dirac comb on Z is periodic, which means that 
there is an integer p > such that Wn+p = Wn for 
all 77 G Z, the diffraction measure is pure point, and 
supported on the lattice Z/p. It is again periodic, at 
least with period 1, but not necessarily with any smaller 
period. As an example, consider the alternating Dirac 
comb with Wn = (—1)"- In this case, 77(777,) = (—1)™ for 
777 G Z, so the autocorrelation is 7 = S27. — 621+1- By 
Poisson's summation formula and elementary properties 
of the Fourier transform (such as the behaviour under 
scaling and the convolution theorem), we have 



cos(27rfc) 



Z/2, 



which leads to the diffraction measure 
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1 — cos(27rfc) 



6z/2 — 6z+l/2 



In this case, the diffraction spectrum is again pure point, 
and consists of Brag peaks of unit intensity at positions 
77, + 1/2, hence on a subset of Z/2. The fundamental 
period of 7 is nevertheless still 1. 

To obtain absolutely or singular continuous compo- 
nents, in line with the classification of Ref. jT^, we thus 
have to go beyond the periodic situation. 

b. Rudin-Shapiro versus Bernoulli. Let us start 
with a deterministic system without periodicity, based 
on the well-known binary Rudin-Shapiro chain. We con- 
sider the corresponding Dirac comb 
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{±1} is defined by the recursion 
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w{n), 

(-l)"+''7«(r7). 



for £ G {0, 1}, 
for £ G {2,3}, 



(2) 



together with the two initial conditions w{0) 1 and 
w{—l) = —1. The resulting system is an aperiodic se- 
quence in 1 and —1, both appearing equally frequent; 
see Fig. [T] for a graphical representation. It has many 
nice properties, such as strict ergodicity and linear patch 
counting complexity, see Ref. ilS and references therein 
for details. In particular, these properties imply that this 
sequence has topological (and metric) entropy 0. 

Since Wj^g is a Dirac comb on Z, the autocorrelation is 
of the form ([T]), with coefficients 
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It follows from the unique ergodicity of the RS sequence 
that all these coefficients (and hence 7) exist, by an ap- 
plication of the ergodic theorem. They are given by 
'7rs(™) ~ ^mO' which is astonishing. The construction 
of a deterministic sequence with vanishing two-point cor- 
relations was the original (and independent) motivation 
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FIG. 1: Central part of the Rudin-Shapiro chain, with open (fuU) circles representing scattering strengths w{n) = 1 {w(n) = — 1), 
respectively. The location of the origin (n — 0) is indicated by the vertical line. 



of Rudin and Shapiro, thus answering a question in the 
theory of Fourier series. 

Let us prove this property by a simple, explicit argu- 
ment. Consider am — Vnsi''^) together with 
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(— 1)" w{n) w(n + m) 



which also exist (by another application of the ergodic 
theorem). Clearly, one has Cq = 1 and = 0, because 
w(ri)^ — 1 for all n S Z. Then, considering to modulo 
4, and splitting the sums in the definition of a„i and hm 
accordingly, the recursion relations ([2]) imply 
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_ l + (-l)'" (-1)'" 7 ill, 

a4m+3 — 4 Clm+l 4 0„j + ^ Orn+l- 

Similarly, one finds 



^4m+l 



0, fo4m+2 

l-(-l)" 



(-ir 



(-1)" 



04m+3 — 4 — 



X + \ bm+1, 



Using the initial data, these recursion relations imply 
that = 5,„ = for all integers m ^ 0. This result 
shows that the autocorrelation and diffraction measures 
of the binary Rudin-Shapiro Dirac comb Wj^g are simply 

7rs = and ^ \, 

where A denotes Lebesgue measure. In other words, the 
diffraction measure is purely absolutely continuous, and 
consists of a constant background (of height 1) only. The 
extinction of all Bragg peaks is due to the balanced choice 
of weights. This is convenient for the theoretical argu- 
ment, but also relevant in practice^^ when disregarding 
thermal displacement. Note that all arguments can be 
extended to mixed spectra. 

Perhaps the most elementary stochastic system is 
based upon the classic coin-tossing (or Bernoulli) exper- 
iment. We consider a stochastic Dirac comb^S on Z, 

nez 

where {Wn)^^^ is a family of independent identically 
distributed (i.i.d.) random variables, with probabilities 
F{Wn = 1) = p and ¥{Wn = -1) = 1 - p, where 
< p < 1. The corresponding (metric) entropy is 



Hip) = -plog(p) - (1 -_p)log(l ~_p), 



(3) 



which satisfies < H{p) < log(2). It attains the extremal 
values for p — and p = 1, where H — is minimal, and 
for p = ^, where H — log(2) is maximal. 

The autocorrelation 73 = X^ngz '7b('^) once 
again a pure point measure that is supported on Z, with 
autocorrelation coefficients given by 



lim 
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for TO G Z. These coefficients almost surely exist, by 
an application of the strong law of large numbers (see 
below), for all m e Z, and satisfy 



?7B(m) = 



m = 



\(2p-l)2, m^O 



This statement can be proved as follows. It obviously 
holds for TO = 0, so consider some fixed to 7^ 0. The prod- 
ucts Zn '■— WnWn+m form a family (Zn)^^^ of identi- 
cally distributed random variables, which take the values 
1 and —1 with probabilities p^ + (f — p)^ and 2p{l—p), re- 
spectively. These new random variables are not indepen- 
dent, but we can split the sum in Eq. ^ into two sums 
(for instance according to even and odd values of [^], the 
largest integer smaller than or equal to — ). The result- 
ing two sums each comprise pairwise independent ran- 
dom variables. An application of the strong law of large 
numbers in its formulation by Etemadi'^^ then shows that 
each sum almost surely converges (as N — > 00) to the ex- 
pectation value of any of the single random variables in 
the sum, which is ^{2p— 1)^. Hence, the diffraction mea- 
sure of the stochastic Dirac comb Wg almost surely is 



i2p-lfS^ + Ap{l-p)X, 



where A again denotes Lebesgue measure. In other 
words, the diffraction spectrum comprises a constant 
background of intensity 4p(l — p) for any value of the 
wave number k and Bragg peaks of intensity {2p — 1)^ 
at integer k. Note that, for the perfectly ordered cases 
p = and p = 1, the background vanishes, while the 
Bragg peaks vanish for the maximally disordered case 
p — ^. At this value of p, the diffraction measure coin- 
cides with that of the Rudin-Shapiro chain. 

This establishes the homometry of the deterministic bi- 
nary Rudin-Shapiro chain (with entropy 0) and the com- 
pletely random Bernoulli chain with p (with entropy 
log(2)), as originally observed in Ref. ,10. Coupling the 
two systems in a suitable way, we now extend this to an 
entire family that covers the intermediate entropy range. 
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c. Bernoullisation. The Bernoulli chain discussed 
above is an example of a completely random and 
interaction-free system. In view of real world examples, 
it is interesting to explore what happens if one imposes 
the influence of coin tossing on the order of a determinis- 
tic system. This can be realised in many different ways. 
Here, we focus on binary sequences and modify them by 
an i.i.d. family of Bernoulli variables. 

Consider a bi-infinite binary sequence S G {±1}^ 
which we assume to be uniquely ergodic. Then, the cor- 
responding Dirac comb ujg = X^nGZ possesses the 
unique (natural) autocorrelation 7^ — X^mGZ ^s('^) '^m 
with the autocorrelation coefficients r]g{m), where 
775(0) = 1 by construction. 

Let (W^„)„gz t>e an i.i.d. family of random variables 
that each take values -1-1 and —1 with probabilities p 
and I — p. The 'Bernoullisation' of LUg is the random 
Dirac comb 



SnWn 5n, 



(5) 



which emerges from by independently changing the 
sign of each Sn with probability 1 — p. Setting Z„ := 
SnWn defines a new family of independent (though in 
general not identically distributed) random variables, 
with values in {±1}. Despite this modification, the auto- 
correlation 7 of w almost surely exists and can be deter- 
mined via its autocorrelation coefficients 77(771) as follows. 
Since one always has 77(0) = 77^(0) = 1, let to 7^ and 
consider, for large N, the sum 
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2A^+ 1 
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WnW, 
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which is split according to the value of {Sn, Sn+m)- Each 
of the four sums can then be handled in the same way 
as in the argument for the Bernoulli chain above, thus 
contributing {2p — 1)^ times the frequency of the corre- 
sponding sign pair. Observing that the overall signs are 
the products SnSn+m, it is clear that, as — > 00, one 
(almost surely) obtains 

r]{m) {2p- ifTJsim) 

for all 771 7^ 0. This shows that the new autocorrelation 
almost surely is 

7= (2p-l)2 7s + 4p(l-p),5o 

where 7^ is the unique autocorrelation of Ug. 

Let us apply this Bernoullisation procedure to the 
Rudin-Shapiro chain. Denote by co the random Dirac 
comb obtained from the Bernoullisation (with parameter 
p) of the binary Rudin-Shapiro chain. Then, the autocor- 
relation measure almost surely exists and reads 7 = ^q. 



independently of p. This means that the random Dirac 
combs oj, even for different values of p, are almost surely 
homometric, and share the purely absolutely continuous 
diffraction measure 7 = A. 

Note that this example explores the fvill entropy range 
of Eq. ([3]): the Bernoulli case (with p — has entropy 
log(2), the maximal value for a binary system, while 
Rudin-Shapiro has entropy 0, and the parameter p in- 
terpolates continuously between the two limiting cases. 
The solution of the corresponding inverse problem is thus 
highly degenerate. Unless additional information is avail- 
able, for instance via higher order correlations, one possi- 
ble strategy could employ a maximum entropy method^, 
singling out the Bernoulli comb. 

Both the Bernoullisation procedure and the spe- 
cific one-dimensional examples immediately generalise to 
higher dimensions by taking direct product structures. In 
particular, the product of d Rudin-Shapiro chains results 
in a deterministic system in d-space with the same purely 
absolutely continuous diffraction measure as the corre- 
sponding coin-tossing model. Consequently, our above 
conclusions extend to this case. This means that one can 
also produce examples with lower rank entropy, which is 
a new phenomenon in dimensions d > 2. 

d. Concluding remarks. Diffraction methods pro- 
vide the most important approach to structure deter- 
mination. The presence of Bragg diffraction clearly in- 
dicates an ordered structure, though the discovery of 
quasicrystals^'^ in the 1980s has shown that pure point 
diffraction occurs in more general systems than just con- 
ventional crystals. To date, the precise atomic structure 
of quasicrystalline alloys is still not completely under- 
stood; there is evidence that entropy plays an important 
role in stabilising quasicrystalline structures, and that 
some disorder may be an inherent feature of these al- 
loys. Like for many ordinary crystals, diffuse scattering 
is present in experimental diffraction patterns of even the 
best known quasicrystals, and there is an increasing effort 
to explore the information contained the diffuse diffrac- 
tion intensity^'^; see Ref. 4 for a recent example. It is 
tempting to draw conclusions about the degree of order 
in a structure on the basis of the observed diffuse scatter- 
ing intensity. However, as our explicit example demon- 
strates, such conclusions have to be considered carefully, 
since the relation between diffuse scattering and disorder 
is far from simple. 
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